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ABSTRACT 


Prototype  operational  FORTRAN  software  is  developed  implementing  a  two-scale 
microwave  sea  surface  scatterometry  model.  A  program  suitable  for  specific 
numerical  testing,  and  another  program  illustrating  its  potential  operational 
utility  in  generating  graphical  visual  aids,  are  also  documented.  Limitations 
of  the  selected  scatterometry  model  are  discussed,  and  suggestions  on  the 
direction  of  future  development  efforts  are  offered. 
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INCORPORATION  OF 
RADAR  SEA  CLUTTER  PREDICTION 
INTO  OPERATIONAL  NAVY  ENVIRONMENTAL  SUPPORT  PRODUCTS: 
PROTOTYPE  SOFTWARE  DEVELOPMENT 


1 .  INTRODUCTION 

It  is  well  known  that  the  normalized  radar  cross  section1  of  the  sea 
surface  displays  great  variability  under  different  environmental  conditions. 
Data  published  by  Guinard,  et  al  (1971),  for  example,  clearly  illustrates  this 
variability;  except  for  a  single  data  set,  they  also  show  a  fairly  regular 
variation  with  respect  to  a  single  parameter- - local  wind  speed.  In  Cheung's 
recent  report  (1988)  and  paper  (1989),  which  review  the  extensive  literature 
relating  to  estimation  of  radar  sea  clutter  based  on  available  (known,  esti¬ 
mated,  or  predicted)  environmental  information,  specific  reference  is  made  to 
a  number  of  theoretical  models  that  have  been  published- -some  of  which  have 
been  reduced  to  computer  code.  Cheung's  rport  gives  reason  to  expect  that 
operational  support  software  appropriate  t  ■■  r  execution  on  small  computers  with 
limited  storage  capacity  and  processing  speed,  could  be  developed  to  predict 
radar  sea  clutter  with  sufficient  fidelity  to  be  tactically  useful. 

In  order  to  assess  the  technical  feasibility  and  potential  utility  of  an 
organic  radar  sea  clutter  estimation  capability,  one  or  more  of  the  available 
theoretical  models  must  be  implemented  in  computer  code  designed  to  facilitate 
embedding  in  prototype,  developmental,  or  even  operational  computer-based 
fleet  support  systems.  This  study  reports  the  successful  coding  of  a  FORTRAN 
subroutine,  termed  "NRCS",  which  implements  a  two-scale  microwave  scattero- 
metry  model  proposed  by  Plant  (1986).  Plant's  model  is  based  on  the  so-called 
"composite  roughness"  theory  that  is  widely  regarded  as  appropriate  for  angles 
of  incidence  in  the  range  bounded  approximately  by  30  and  70  degrees  ("inter¬ 
mediate"  angles).  The  subroutine  performs  neither  READ  nor  WRITE  operations 
that  would  preclude  its  integration  into  specialized  systems  that  lack  support 
for  ordinary  FORTRAN  input/output  units,  so  a  test  program  "TNRCS"  is  provided 
to  facilitate  testing  for  user-specified  numerical  inputs.  Program  "DNRCS"  is 
also  provided  to  demonstrate  a  concept  of  how  the  subroutine  could  be  employed 
to  generate  a  tactically  useful  computer-generated  visual  aid. 


2.  CODING  AND  TESTING  AN  EMBEDDED  NORMALIZED  RADAR  CROSS  SECTION  SUBROUTINE 

The  author  was  fortunate  to  be  able  to  start  with  the  research-oriented 
FORTRAN  program  SIGMOD  originally  developed  by  Plant  and  subsequently  modified 
by  Plant  and  Cheung.  That  program  successfully  implements  the  theory  devel- 


1The  normalized  radar  cross  section  a0  of  a  region  of  the  ocean  surface  is  the 
ratio  of  radar  cross  section  to  its  physical  area.  Thus,  if  the  radar  return 
from  a  10  m2  region  were  equal  to  that  from  a  perfect  reflector  1  m2  in  area, 
the  normalized  radar  cross  section  would  equal  0.1,  or  -10  dB  in  logarithmic 
terms  (dB  -  10  log10  a0  )  . 
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oped  in  Plant  (1986),  but  is  unsuitable  for  embedding  into  operational  soft¬ 
ware;  it  is,  itself,  a  main  program  that  accepts  input  data  from  a  keyboard, 
generates  an  array  of  normalized  radar  cross-section  values  over  a  range  of 
parameter  values,  outputs  summary  information  to  the  monitor,  and  writes 
tabular  information  to  two  different  disk  files. 

Subroutine  NRCS  is  designed  to  return  a  single  value  of  normalized  radar 
cross  section  for  a  given  set  of  parameters  describing  environmental  condi¬ 
tions  and  microwave  radiation  characteristics.  Because  it  performs  no  input/ 
output  operations  of  its  own,  it  can  be  embedded  in  any  system  capable  of 
executing  FORTRAN  subroutine  calls  and  its  execution  is  controlled  by  the 
invoking  program.  Appendix  A  lists  the  subroutine  code,  which  includes 
comments  defining  the  arguments.  Appendix  B  is  a  listing  of  program  TNRCS, 
which  was  used  to  test  NRCS  output  against  tables  generated  by  program  SIGMOD. 

Because  NRCS  does  not  record  intermediate  quantities  used  in  the 
computation,  SIGMOD  remains  the  preferred  tool  for  research  study  of  the 
behavior  of  Plant's  model  as  a  function  of  various  inputs. 


3.  ILLUSTRATING  THE  POTENTIAL  TACTICAL  UTILITY  OF  SEA  CLUTTER  ESTIMATION 

Program  DNRCS  was  written  to  demonstrate  the  potential  tactical  utility 
of  an  organic  capability  to  estimate  radar  sea  clutter.  It  accesses  a  data 
file  containing  wind  fields  over  a  region  of  the  ocean  obtained  from  an  opera¬ 
tional  weather  product  from  the  Fleet  Numerical  Oceanography  Center.  Calls  to 
NRCS  are  executed  for  each  of  the  grid  points,  and  three  graphical  displays 
are  then  generated.  Figures  1  and  2  show,  respectively,  absolute  and  loga¬ 
rithmic  contour  plots  from  DNRCS  of  normalized  radar  cross  section.  Of 
greater  interest  would  be  something  like  Figure  3,  which  displays  estimated 
signal- to-noise-ratio  of  a  target  of  interest  with  respect  to  a  background  of 
sea  clutter.  Appendix  C  is  a  program  listing  of  DNRCS. 


4.  ENGINEERING  A  FLEXIBLE  SOFTWARE  PACKAGE  TO  FACILITATE  EMBEDDING  OF  SEA 

CLUTTER  PREDICTION  IN  OPERATIONAL  SYSTEMS 

One  can  imagine  a  number  of  operational  scenarios  in  which  radar-depen¬ 
dent  surveillance  sensors  and  weapons  could  be  affected  by  sea  clutter;  hence, 
the  ability  to  estimate  current  effects  of  clutter,  and  to  predict  them,  could 
be  advantageous  to  a  tactical  commander.  Appendix  D  sketches  out  a  simple 
conceptual  software  package,  Environmental_Sea_Clutter ,  in  pseudocode 
mimicking  the  Ada  computer  language,  which  would  provide  three  different 
levels  of  information  on  environmentally- related  sea  clutter.  The  lowest- 
level  procedure,  NRCS,  would  estimate  normalized  radar  cross  section  as  a 
function  of  environmental  parameters  that  could  reasonably  be  expected  to  be 
available  aboard  a  naval  platform  either  from  direct  observation  or  from 
routine  external  environmental  support  sources.  The  FORTRAN  subroutine  of 
Appendix  A  can  be  considered  a  first  prototype  of  this  procedure,  but  a  model 
of  somewhat  higher  fidelity  would  probably  be  needed.  Treatment  of  near¬ 
normal  and  near-grazing  incidence  angles  will  also  have  to  be  included. 
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Figure  1.  Contour  plot  of  normalized  radar  cross  section  o0  generated 
by  program  DNRCS . 


It  is  unlikely  that  the  output  of  NRCS  would  be  of  direct  interest  to 
operational  commanders,  whose  concern  would  inevitably  focus  on  the  capability 
of  a  given  radar-based  system  (his  own  or  an  adversary's)  to  perform  effec¬ 
tively.  A  higher- level  procedure,  ERCS ,  would  be  required  to  estimate  envi¬ 
ronmental  radar  cross  section  based  on  NRCS  and  parameters  characterizing  the 
radar  itself  (e.g.  pulse  duration,  beamwidth)  and  its  location  with  respect  to 
a  target  of  interest  (e.g.  altitude,  azimuth).  Procedure  TSNR  would  then  be 
employed  to  estimate  signal-to-noise  ratio  of  the  target  with  respect  to 
environmentally- generated  clutter,  based  on  the  ERCS  and  the  target's  own 
radar  cross  section.  The  latter  procedure  could  be  used  to  generate  displays 
along  the  lines  of  Figure  3,  and  could  also  serve  as  the  basis  for  higher- 
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Figure  2.  Contour  plot  of  normalized  radar  cross  section  in  dB 
(dB  -  10  log10  trQ)  .  Same  field  as  Figure  1. 

level  displays  of  where  certain  targets  would  or  would  not  be  detectable  by 
the  radar  in  question. 

The  software  engineering  effort  to  develop  Environmental_Sea_Clutter 
would  probably  identify  other  modules  for  that  generic  "toolbox"  package. 

5.  RECOMMENDATIONS  FOR  FUTURE  STUDY 

Efforts  to  continue  the  line  of  development  initiated  in  this  study  can 
logically  be  divided  into  three  categories:  (1)  identification  of  where,  and 
to  what  extent,  a  Navy  need  exists  for  an  "rganic  radar  sea  clutter  ^uimation 


Figure  3.  Estimated  signal - to-noise  ratio  (SNR)  in  dB  of  a  target  of 
interest  with  respect  to  a  background  of  sea  clutter.  Same  field  as 
Figure  1. 

capability- -  to  include  conceptual  development  of  how  this  kind  of  information 
could  be  displayed,  or  otherwise  used,  in  existing/developmental/proposed 
computer-based  systems;  (2)  analysis  to  identify  what  requirements  the  above 
needs -oriented  study  imposes  on  a  "toolbox"  package  along  the  lines  of 
Appendix  D;  and  (3)  development  of  the  lowest-level  module  "NRCS",  with 
fidelity  sufficient  to  satisfy  operational  requirements  without  violating  the 
constraints  that  only  readily  available  data  can  be  used  and  that 
computational  load  not  exceed  the  organic  resources  of  operational  platforms. 

The  first  and  second  categories  of  effort  together  constitute  a  top-down 
exploratory  development  (6.2)  program  that,  if  successful,  would  lead  to 
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transitions  into  advanced  development  (6.3)  programs  and/or  upgrades  to 
existing  fleet  systems.  Hopefully,  the  mainline  6.2  effort  would  find  support 
us  part  of  one  of  the  Office  of  Naval  Technology’s  exploratory  development 
block  programs.  Possibly,  a  case  could  be  made  now  to  develop  an  application 
for  some  particular  developmental  or  operational  system,  in  which  case  6.3  (or 
higher  category)  funding  should  be  sought  from  the  appropriate  sponsor. 

In  the  case  of  the  third  category- -  the  effort  to  implement  the  required 
underlying  mathematical  capability- - there  is  probably  a  need  for  a  top-down 
effort  to  define  the  required  capabilities  as  well  as  a  bottom-up  effort  to 
determine  whether  adequate  theoretical  models  are  available.  Requirements 
definition  and  selection  of  models  to  be  implemented  would  seem  to  belong  in 
the  6.2  arena,  whereas  6.1  funding  should  be  sought  for  whatever  fundamental 
research  is  needed  to  address  deficiencies  in  the  current  "state-of-the-art”. 
Since  it  appears  highly  unlikely  that  a  single  model  will  be  valid  over  the 
entire  possible  range  of  incidence  angles  (0  to  90  degrees) ,  the  module  that 
finally  emerges  from  this  effort  will  need  to  exhibit  smooth  behavior  in  the 
transition  regions  between  models  appropriate  to  near-normal,  intermediate, 
and  grazing  incidence  angles. 

There  is  a  definite  liklihood  that  additional  fundamental  research  will 
be  needed.  Plant,  1986,  states  that  his  proposed  model  is  incomplete  and  is 
mainly  intended  to  indicate  which  physical  variables  require  additional 
investigation  for  algorithm  development.  The  principal  deficiency  of  this 
model  is  probably  the  fact  that  the  long-wave  ("swell")  component  of  the  local 
surface  wave  motion  is  assumed  to  have  reached  equilibrium  for  the  specified 
wind  parameters.  (Note  that  Subroutine  NRCS  of  Appendix  A  does  not  allow  for 
specification  of  the  direction  of  long-wave  motion.)  Plant  addresses  this  and 
other  deficiencies  in  the  conclusion  of  his  paper.  It  is  worth  noting  that 
Plant's  work  is  focused  on  the  problem  of  estimating  the  values  of  certain 
atmospheric  and  oceanic  parameters  based  on  satellite-mounted  scatterometry 
data,  which  is  the  inverse  of  our  desire  to  estimate  sea  clutter  as  a  function 
of  environmental  parameters;  perhaps  a  research  effort  explicitly  pursued  from 
the  latter  perspective  would  be  fruitful. 
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APPENDIX  A.  NRCS  LISTING  PAGE  1 

-  03-16-90 

16:40:31 

Line/  Source  Line  Microsoft  FORTRAN  Optimizing  Compiler  Version  4.01 

1  C  SUBROUTINE  NRCS  VERS  16  MAR  90 

2  C 

3  C  MODULE  TO  COMPUTE  THE  NORMALIZED  RADAR  CROSS  SECTION  OF  THE  OCEAN 

4  C  SURFACE  AS  A  FUNCTION  OF  ENVIRONMENTAL  FACTORS 

5  C 

6  C  ARGUMENTS:  [ (R)  DENOTES  REAL,  (I)  DENOTES  INTEGER) 

7  C  INPUTS: 


8 

C 

WNM 

(R) 

RADAR  WAVENUMBER  (INVERSE  CM) 

9 

C 

RANGE:  WNM  >  0 

10 

c 

POL 

(I) 

POLARIZATION  (DIMENSIONLESS  INDEX) 

11 

c 

POL  -  0  IMPLIES  VERTICAL,  POL  -  1 

IMPLIES  HORIZONTAL 

12 

c 

PHI 

(R) 

ANGLE  OF  INCIDENCE  (RADIANS) 

13 

c 

RANGE:  0  <  PHI  <  PI/2 

14 

c 

CHI 

(R) 

AZIMUTH  WITH  RESPECT  TO  WIND  (RADIANS) 

15 

c 

RANGE:  0  <-  CHI  <  2*PI 

16 

c 

SLOPE 

(R) 

MEAN  LONG-WAVE  SLOPE 

17 

c 

RANGE:  0  <  SLOPE  (OR  SLOPE  »  -1.0 

FOR  COMPUTED  DEFAULT 

18 

c 

WIND 

(R) 

WIND  SPEED  (M/SEC) 

19 

c 

RANGE:  0  <*  WIND  <»  26  M/SEC 

20 

c 

OUTPUTS : 

21 

c 

SIG0 

(R) 

NORMALIZED  RADAR  CROSS  SECTION  (DIMENSIONLESS) 

22 

c 

RANGE:  0  <  SIGO 

23 

c 

SIGDB 

(R) 

•  ••(?)••• 

24 

c 

CA0DB 

(R) 

•  »•(?)••• 

25 

c 

IRET 

(I) 

RETURN  CODE: 

26 

c 

0  - 

SUCCESSFUL  RETURN 

27 

c 

-1  - 

WNM 

OUT  OF  RANGE 

28 

c 

-2  - 

POL 

VALUE  ILLEGAL 

29 

c 

-3  - 

PHI 

OUT  OF  RANGE 

30 

c 

-4  - 

CHI 

OUT  OF  RANGE 

31 

c 

-5  - 

A 

OUT  OF  RANGE 

32 

c 

-6  - 

SLOPE  OUT  OF  RANGE 

33 

c 

-7  - 

WIND  OUT  OF  RANGE 

34 

c 

1  - 

UW/C0  >1.7  (INTERMEDIATE  COMPUTATION 

FAILURE  CONDITION) 

35 

c 

99  - 

UNSPECIFIED  FAILURE 

36 

c 

37 

c 

EMPIRICAL  PARAMETERS: 

38 

c 

A 

<R) 

GROWTH  RATE  PARAMETER  OF  PAGE  10,737  (DIMENSIONLESS) 

39 

c 

EPS 

(R) 

CROSSWIND/UPWIND  SLOPE  RATIO  (DIMENSIONLESS) 

40 

c 

AA 

(R) 

PARAMETER  "AO"  OF  EQ  (49) 

41 

c 

BB 

(R) 

PARAMETER  "Al"  OF  EQ  (49) 

42 

c 

CC 

(R) 

PARAMETER  "A2"  OF  EQ  (49) 

43 

c 

44 

SUBROUTINE  NRCS (WNM, POL , PHI , CHI , SLOPE , WIND , 

45 

1 

IRET, SIGO, SIGDB, CAODB) 

46 

c 

47 

c 

ARITHMETIC  CONSTANTS: 

48 

REAL 

PI , HALFPI , THPI , TWOPI 

49 

PARAMETER 

(  PI  -  3.1415927,  HP I  -  PI/2.0, 

50 

1 

THPI  -  3 . 0*HPI ,  TWOPI  -  2 . 0*PI  ) 

51 

c 

52 

c 

EMPIRICAL  CONSTANTS  (PER  PAGE  10,741): 

53 

REAL  , 

A,  AA, 

BB,  CC 
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Line#  Source  Line 


PAGE  2 
03-16-90 
15:40:31 


54 

55  C 

56  C 

57 

58 

59 

60 

61  C 

62 

63 

64  C 

65  C 

66 

67 

68 

69 

70  C 

71  C 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 


Microsoft  FORTRAN  Optimizing  Compiler  Version  4.01 
PARAMETER  (  A  -  0.05,  EPS  «  0.4,  AA  -  1.0,  BB  -  0.05,  CC  -  -1.0  ) 

TYPE  DECLARATIONS: 

INTEGER  POL 

REAL  WNM, PHI, CHI, SLOPE, WIND 

INTEGER  IRET 

REAL  SIG0, SIGDB, CA0DB 

REAL  SSIN , CCOS , TTAN , CCOT , CD , WNV , CWIND , LWSLP 
REAL  PMV , PM2 V , PMH , PM2H , G , GV , GH , T3 

INITIALIZATION: 

IRET  «  99 
SIG0  «  -1.0 
SIGDB  «  1000000.0 
CA0DB  »  1000000.0 

ENSURE  INPUTS  ARE  WITHIN  ALLOWED  RANGE: 

IF  (WNM. LE. 0.0)  THEN 
IRET  -  -1 
GOTO  900 
ENDIF 

IF  (  (POL.NE.O) .AND. (POL.NE.l)  )  THEN 
IRET  -  -2 
GOTO  900 
ENDIF 

IF  (  (PHI.LE.0.0) .OR. (PHI.GE.HPI)  )  THEN 
IRET  -  -3 
GOTO  900 
ENDIF 

IF  (  (CHI.GE.0.0) .AND. (CHI.LT.TWOPI)  )  THEN 
IF  (  (CHI.LT.HPI) .OR. (CHI.GT.THPI)  )  THEN 
THETA  ■  CHI 
ELSE 

THETA  -  CHI  -  PI 
ENDIF 
ELSE 

IRET  -  -4 
GOTO  900 
ENDIF 

IF  (A.LE.0.0)  THEN 
IRET  -  -5 
GOTO  900 
ENDIF 

IF  (  (SLOPE. LE. 0.0) .AND. (SLOPE. NE. -1.0)  )  THEN 
IRET  -  -6 
GOTO  900 
ENDIF 

IF  (  (WIND. LT. 0.0) .OR. (WIND. GE. 26.0)  )  THEN 
IRET  ■  -7 
GOTO  900 
ENDIF 

SSIN  -  SIN(PHI) 
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Line# 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 
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CCOS  -  COS (PHI) 

TTAN  »  SS IN /CCOS 
CCOT  -  1/TTAN 
CWIND  -  100 . 0*WIND 
C 

C++++++++++++++++++COMPUTE  WIND  AND  WAVE  PARAMETERS+++++++++++++++++++++ 
C 

C  DRAG  COEFFICIENT  "CD"  -  EQ  ;46): 

IF  (WIND. LE. (4.0) )  THEN 

CD  «  0.00114 

ELSEIF  ( WIND. LT. 10.0)  THEN 
CD  -  0.00114 

ELSEIF  (WIND. LT. (26.0))  THEN 
CD  -  0.00049  +  0 . 000065*WIND 
END  IF 
C 

C  FRICTION  VELOCITY  "DSTAR": 

USTAR  -  CWIND "SQRT ( CD ) 

C 

WNW-  2 . 0*WNM*SSIN 

CO  «  SQRT(  9B1/WNW  +  74.0*WNW  ) 

C 

C  "UW"  -  EQ  (36): 

UW  ■  0. 6 "USTAR  -  0 . 084 "USTAR* ALOG (  1.0  +  92.0/WNW  ) 

IF  (UW.GT. (1.7*C0) )  THEN 
IRET  -  1 
GOTO  900 
ENDIF 
C 

C  "PM0"  -  EQ  (49): 

PM0  -  AA  +  BB*USTAR  +  CC*WNW 
C 

C  -  CO  +  UW*COS (THETA) 

AFR  -  WNW*C 

GTK  -  981.0  +  222.0*WNW*WNW 
CG  -  GTK/ ( 2 . 0*WNW*C0 ) 

1  +  (UW+7 . 728*DSTAR/  (WNW+92 . 0) ) *COS (THETA) 

DCGDK  -  —GTK* *2/ (4.0*WNW**3*C0**3)  +  222.0/C0 
1  +  7.728*OSTAR*COS  (THETA)  *(92. /(WNW*  (WNW+92. )))/ (WNW+92. ) 

CGC  -  CG/C 

BET0  -  0. 04*USTAR**2*AFR/C**2 
C 

C  "BETW"  -  EQ  (38): 

BETW  »  0.04*USTAR**2/(C0+0.5*UW)**2 
C 

C  "R"  -  EQ  (C4) : 

R  -  UW/ (C0+0. 5*UW) 

C 

C  EXPANSION  COEFFICIENTS  IN  "R"  -  EQ  (C8): 

R1  »  0.64  -  0. 36*R  +  0.25*R**2 

R2  -  0.42  -  0. 58*R  +  0.34*R**2 

R4  »  0.08  +  0. 08*R  -  0.16*R**2 

R6  ■  0.04  ♦  0.04*R  +  0.02*R**2 
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C++++++++++++++++++++++COMPUTE  mp  &  Ti  VALUES++++++++++++++++++++-t.+-*.++++ 
C 

GH  -  (CCOS/(0.111*CCOS+1.0))**4 

GV  -  CC0S**4* ( 1 .0+SSIN**2 ) **2/ (CCOS+O. Ill ) **4 

IF  (POL.EQ.0)  THEN 

PMV  -  2.0*SIN(2*PHI)/(1.0+SSIN**2)  -  4.0*TTAN 
1  +  4.0*SSIN/ (CCOS+O.  ill) 

PM2V  -  PMV* *2  -  4.0*(3.0*CCOS**2-2.0) / ( 1+SSIN**2 ) **2 

1  +  4.0*  ( 1.0+0. lll*CCOS )  /  (CCOS+O .  Ill )  **2 

2  -  4.0/CCOS**2 

PM2P  -  PM2V 

PMP  -  PMV 
G  «  GV 

T3  -  2 . 0* ( SQRT ( GH/GV)  -  0.5)/SSIN**2 
ELSEIF  (POL.EQ.l)  THEN 

PMH  -  -4.0*TTAN  +  0.444*SSIN/{0.111*CCOS+1.0) 

PM2H  -  PMH**2  -  4.0/CCOS**2 
1  +  0. 444* (0. 111+CCOS) / (0. lll*CCOS+l. 0) **2 

PM2P  «  PM2H 
PMP  -  PMH 
G  -  GH 

T3  -  2.0* (SQRT (GV/GH)  -  0.5)/SSIN**2 
ELSE 

IRET  -  99 
GOTO  900 
ENDIF 

T2  «  PM2P/2.0  -  PMP*CCOT* ( 2 . 0*CGC+2 . 0 ) 

1  +  ( 4 . 0*CCOT**2+1 . 0 ) *CGC  +  3.0*CCOT**2  +  1.0 

2  +  (3.0*CGC**2  -  WNW*DCGDK/C) *CCOT**2 
Tl  -  PMP  +  CCOT*  (WNW*CC/PM0-2 .0*CGC-2 .0 ) 

C 

C++++++++++++COMPUTE  SLOPE  AND  COEFFICIENTS  A1++++++++++++++++++++++++++ 
C 

IF  (SLOPE. EQ. -1.0)  THEN 

LWSLP  ■  0.008  +  0. 0000156 *CWIND 
ELSE 

LWSLP  -  SLOPE 
ENDIF 

A00  -  BETW*R1 

A01  ■  LWSLP* (  0 . 5*BETW*R1* (T2+T3 ) 

1  +  0. 25*BETW*R2* (1.0-2 . 0*EPS ) * (T2-T3 ) 

2  -  0. 5*CGC*PM0/AFR  +  0. 5*WNW*CC/AFR  ) 

A0  -  A00  +  A01 

A1  -  LWSLP* (  0. 5*PM0*Tl*BETW*(Rl+0. 5*R2 )  ) 

A20  -  BETW*R 2 

A21  -  LWSLP*  (0. 5*BETW*R2* (T2+T3 ) 

1  +  0 . 2S*BETW* ( 2 .0*R1-R4) • (1.0-2.0*EPS)* (T2-T3 ) 

2  -  (2 . 0+0 . 5*CGC) *PM0/AFR  +  0. 5*WNW*CC/AFR) 

A2  -  A20  +  A21 

A2P  -  -LWSLP*  (PM0*UW*  SIN  (THETA)  *0. 5/  (C*AFR) ) 

A3  -  LWSLP*  (0. 25*PM0*T1*BETW*  (R2-R4 ) ) 

A40  -  -BETW*R4 


A00  - 
A01  - 
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213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 


A41  -  LWSLP* (  -0. 5*BETW*R4* (T2+T3) 

1  +  0 . 25*BETW* (R2+R6) *(1.0-2. 0*EPS ) * (T2-T3 )  ) 

A4  ■  A40  +  A41 

C 

C+++++++++++++++++EVALUATE  SIGMAO  AND  TAKE  LOG++++++++++++++++++++++++++ 
C 

COEF  -  PI*A*G/SSIN**4 

SIG0  «  COEF* (  A0  +  Al*COS ( CHI ) 

1  +  A2  *COS (2 . 0*CHI )  +  A3*COS<3.0*CHI) 

2  +  A2P*SIN( 2 . 0*CHI )  +  A4*COS<4.0*CHI)  ) 

RATIO  -  (  A01+A1*COS (CHI )  +  A21*COS(2.0*CHI)  +  A3*COS(3.0*CHI) 

1  +  A41*COS(4.0*CHI)  +  A2P*SIN(2.0*CHI)  ) 

2  /(  A00  +  A20*COS (2 . 0*CHI )  +  A40*COS(4.0*CHI)  ) 

IF  (SIG0.GT.0.0)  THEN 

SIGDB  -  10. 0*ALOG10 (SIGO) 

ELSE 

SIGDB  ■  -99.0 
ENDIF 

CAO  »  COEF*AO 
IF  (CAO. GT. 0.0)  THEN 

CAODB  -  10.0*ALOG10(CA0) 

ELSE 

CAODB  «  -99.0 
ENDIF 
C 

C* ********************************************************************** 
IRET  -  0 
900  CONTINUE 
RETURN 
END 


NRCS  Local  Symbols 


Name 

Class 

Type 

Size 

Offset 

CAODB  . 

0006 

SIGDB  . 

000a 

SIGO . 

OOOe 

IRET . 

0012 

HIND . 

0016 

SLOPE  . 

001a 

CHI . 

OOle 

PHI . 

0022 

POL . 

0026 

WNM . 

002a 

A2P . 

REAL *4 

4 

0000 

R1 . 

REAL*4 

4 

0004 

C . 

REAL* 4 

4 

0008 

R2 . 

REAL* 4 

4 

000c 

T1 . 

REAL* 4 

4 

0010 

T2 . 

REAL *4 

4 

0014 

R4 . 

•  •  •  •  •  local 

REAL* 4 

4 

0018 
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NRCS  Local  Symbols 
Name 


T3.  . 
G  .  . 
CD.  . 
R6.  . 
CG.  . 
BETO. 
CGC  . 
PMO  . 
OSTAR 
GH.  . 
R  .  . 
LWSLP 
PM2H. 
AFR  . 
COEF. 
DCGDK 
GV.  . 
PM2P. 
A00  . 
A01  . 
A20  . 
A21  . 
A40  . 
PMH  . 
PM2V. 
A41  . 
GTK  . 

CCOS. 

CCOT. 
OW.  . 
PMP  . 
AO.  . 
Al.  . 
BETW. 
A2.  . 
CO.  . 
PMV  . 
CAO  . 
A3.  . 
A4.  . 
CWIND 
THETA 
TTAN. 
WNW  . 
SSIN. 
RATIO 


Class 

Type 

Size 

Offset 

local 

REALM 

4 

001c 

local 

REALM 

4 

0020 

local 

REALM 

4 

0024 

local 

REALM 

4 

0028 

local 

REALM 

4 

002c 

local 

REALM 

4 

0030 

local 

REALM 

4 

0034 

local 

REALM 

4 

0038 

local 

REALM 

4 

003c 

local 

REALM 

4 

0040 

local 

REALM 

4 

0044 

local 

REALM 

4 

0048 

local 

REALM 

4 

004c 

local 

REALM 

4 

0050 

local 

REALM 

4 

0054 

local 

REALM 

4 

0058 

local 

REAL* 4 

4 

005c 

local 

REAL *4 

4 

0060 

local 

REALM 

4 

0064 

local 

REALM 

4 

0068 

local 

REALM 

4 

006c 

local 

REALM 

4 

0070 

local 

REALM 

4 

0074 

local 

REALM 

4 

0078 

local 

REALM 

4 

007c 

local 

REAL* 4 

4 

0080 

local 

REALM 

4 

0084 

local 

REALM 

4 

0088 

local 

REALM 

4 

008c 

local 

REALM 

4 

0090 

local 

REALM 

4 

0094 

local 

REALM 

4 

0098 

local 

REAL* 4 

4 

009c 

local 

REALM 

4 

OOaO 

local 

REALM 

4 

00a4 

local 

REALM 

4 

00a8 

local 

REALM 

4 

OOac 

local 

REALM 

4 

OObO 

local 

REALM 

4 

00b4 

local 

REAL* 4 

4 

00b8 

local 

REALM 

4 

OObc 

local 

REALM 

4 

OOcO 

local 

REALM 

4 

00c4 

local 

REALM 

4 

00c8 

local 

REALM 

4 

OOcc 

local 

REAL* 4 

4 

OOdO 
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Symbolic  Constant  Type  Value 

PI . REALM  3 . 1415927E+000 

THPI . REALM  4.7123890E+000 

TWOPI . REALM  6.2831855B+000 

HPI  .  REALM  1 . 5707964E+000 

A . REAL* 4  5 . 0000001E-002 

AA .  REAL*4  1 . OOOOOOOE+OOO 

BB . REAL*4  5 . 0000001E-002 

CC . REALM  -1 .  OOOOOOOE+OOO 

EPS . REALM  4.0000001E-001 


Global  Symbols 


Name 


Class  Type 


Size  Offset 


NRCS 


FSUBRT  *** 


***  0000 


Code  size  «  0d73  (3443) 
Data  size  -  00b4  (180) 
Bss  size  ■  00d4  (212) 

No  errors  detected 
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1  PROGRAM  TNRCS 

2  C 

3  C  ARITHMETIC  CONSTANTS: 

4  REAL  PI , HALFPI , THPI , TWOPI 

5  PARAMETER  (  PI  -  3.1415927,  HPI  -  PI/2.0, 

6  1  THPI  -  3.0*HPI,  TWOPI  ■  2.0*PI  ) 

7  C 

8  REAL  WNM,ANGINC,ANGAZ, SLOPE, WIND 

9  INTEGER  POL 

10  C 

11  WRITE (6, 6000) 

12  6000  FORMAT (IX, '*  *  *  N  R  C  S  Test  Program  *  *  *',/) 

13  C 

14  1  CONTINUE 

15  PRINT  *, 'ENTER  THE  MICROWAVE  NUMBER  IN  CM-1: ' 

16  READ  *,  WNM 

17  C 

18  5  CONTINUE 

19  PRINT  *, 'ENTER  THE  POLARIZATION  (V-0,H»1) : ' 

20  READ  *,  POL 

21  IF  ( POL .  NE .  0 .  AND .  POL .  NE .  1 )  GO  TO  5 

22  C 

23  PRINT  *, 'ENTER  THE  INCIDENCE  ANGLE  (DEG):' 

24  READ*,  ANGINC 

25  C 

26  PRINT  *, 'ENTER  THE  WIND  SPEED  (M/S),  AZIMUTH  ANGLE  (DEG):' 

27  READ  *,  WIND,  ANGAZ 

28  C 

29  PRINT  *, 'Enter  Mean-Square  Long-Wave  Slope  (or  -1  for  default):' 

30  READ  *,  SLOPE 

31  C 

32  PHI  ■  PI* ANGINC/ 180.0 

33  CHI  -  PI* ANGAZ/ 180.0 

34  CALL  NRCS(WNM, POL, PHI, CHI, SLOPE, WIND, 

35  1  IRET , SIG0 , SIGDB , CA0DB ) 

36  C 

37  WRITE  (*,6001)  IRET,  SIGO,  SIGDB,  CAODS 

38  6001  FORMAT (IX,  'IRET  13,'  SIGO  -',  IP,  G12.5,  OP, 

39  1  '  SIGDB  F7.3, 

40  2  '  CAODB  F7 . 3 , / ) 

41  900  CONTINUE 

42  PRINT  *,'  RUN  MODEL  AGAIN?  (TYPE  0  TO  REPEAT,  1  TO  STOP)  ' 

43  READ  *,NN 

44  IF  (NN.EQ.O)  GO  TO  1 

45  END 

main  Local  Symbols 

Name  Class  Type  Size  Offset 

SLOPE . local  REAL*  4  4  0002 

SIGO . local  REAL*  4  4  0006 
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main  Local  Symbols 

Name 

Class 

CHI . 

.  .  local 

NN •*•••••••• 

.  .  local 

.  .  local 

SIGDB  •••«••«• 

.  .  local 

POL  **••••••• 

.  .  local 

ANGINC . 

.  .  local 

ANGAZ  •••••••• 

.  .  local 

WIND ••••••««» 

.  .  local 

WNH  ••••••«•• 

.  .  local 

IRET . 

.  .  local 

CA0DB  . 

.  .  local 

Type 

Size 

Offset 

REALM 

4 

000a 

INTEGERM 

4 

OOOe 

REALM 

4 

0012 

REALM 

4 

0016 

INTEGER* 4 

4 

001a 

REALM 

4 

OOle 

REALM 

4 

0022 

REALM 

4 

0026 

REALM 

4 

002a 

INTEGERM 

4 

002e 

REALM 

4 

0032 

Symbolic  Constant 


Type 


Value 


PI.  . 
THPI . 
TWOPI 
HP  I  . 


REALM 

REALM 

REALM 

REALM 


3 . 1415927E+000 
4.7123890E+000 
6.2831855E+000 
1.S707964E+000 


Global  Symbols 


Name  Class 

. extern 

main . FSUBRT 


Type 

Size 

offset 

*** 

*** 

*** 

*** 

*  *  * 

0000 

Code  size  *  0165  (357) 
Data  size  *  0179  (377) 
Bss  size  »  0036  (54) 


No  errors  detected 
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1 

PROGRAM  DNRCS 

2 

C 

VERS 

16  MAR  1990 

3 

c 

4 

c 

ARITHMETIC  CONSTANTS: 

5 

REAL  PI , HALFPI , THPI , TWOPI 

6 

PARAMETER  (  PI  -  3.1415927,  HPI  -  PI/2.0, 

7 

1  THPI  -  3 . 0*HPI ,  TWOPI  «  2.0*PI  ) 

8 

c 

9 

REAL  WNM ,  ANGINC ,  ANGLE ,  ANGAZ ,  SLOPE ,  RDAREA ,  TGTSIG 

10 

INTEGER  POL 

11 

c 

12 

c 

Z  CONTAINS  THE  VALUES  TO  BE  PLOTTED. 

13 

c 

14 

c 

15 

REAL  Z(21,25),  UWIND (21, 25 ) ,  VWIND(21,25) ,  DZ(21,25), 

SNR(21 , 25 ) 

1  c 

c 

17 

c 

Define  some  default  colors  to  use 

18 

c 

19 

REAL  COLS 

20 

INTEGER  NCOLI ,  COLA,  NPIC,  SLOTS,  COLI 

21 

REAL  RED ( 6 ) ,  BLUE(6),  GREEN ( 6 ) 

22 

DATA  RED/1.,1.,0.,0.,1.,0./ 

23 

DATA  GREEN/1.,0.,0.,1.,1.,1./ 

24 

DATA  BLUE/1.,0.,1.,0.,0.,1./ 

25 

c 

26 

c 

SPECIFY  COORDINATES  FOR  PLOT  TITLES.  ON  AN  ABSTRACT  GRID 

WHERE 

27 

c 

THE  INTEGER  COORDINATES  RANGE  FROM  0.0  TO  1.0,  THE  VALUES 

TX  AND  TY 

28 

c 

DEFINE  THE  CENTER  OF  THE  TITLE  STRING. 

29 

c 

30 

DATA  TX/ . 3955/,  TY/.9765/ 

31 

c 

32 

1  CONTINUE 

33 

PRINT  *, 'Enter  the  Microwave  Number  in  cm-1:' 

34 

READ  *,  WNM 

35 

c 

36 

5  CONTINUE 

37 

PRINT  *, 'Enter  the  Polarization  (V*0,H*1):' 

38 

READ  *,  POL 

39 

IF  ( POL . NE . 0 . AND . POL . NE . 1 )  GO  TO  5 

40 

c 

41 

PRINT  *, 'Enter  the  Incidence  Angle  (deg):' 

42 

READ*,  ANGINC 

43 

c 

44 

PRINT  *, 'Enter  Look  Direction  (deg):' 

45 

READ  *,  ANGLE 

46 

c 

47 

c 

open  GKS  and  active  the  display 

48 

c 

49 

call  xgbeg 

50 

call  xtcopt( 'OPNWKS' ,  1.0) 

51 

c 

52 

c 

53 

c 

Determine  the  number  of  color  slots  we  can  fill  (max  we 
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54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 


C  want  is  5).  Note  that  some  displays  do  not  have  loadable 
C  color  tables  and  the  XTCCOL  will  have  no  effect. 

C 

CALL  XTCOPT ( ' COLAVL ' , COLS ) 

NCOLI  -  COLS 
IF  (NCOLI. EQ.0)  THEN 
C  Monochrome  so  dont  set  any  colors 
SLOTS  «  0 
ELSE 

SLOTS  -  MIN0 (NCOLI-1 , 6) 

CALL  XTCCOL (1,  SLOTS,  RED,  GREEN,  BLUE) 

END  IF 
C 

C  Set  the  text  color  to  index  1. 

C 

CALL  XTCOPT ( 'TCOLOR' ,1.0) 

C 

C  Set  the  size  of  labels 
C 

CALL  CTROPT{ 'ISIZEL' ,2.0) 

CALL  CTROPT( 'ISIZEM' ,3.0) 

C 

C  FILL  TWO  DIMENSIONAL  ARRAY  TO  BE  PLOTTED 
C 

C  GET  WIND  FIELDS: 

OPEN  ( UNIT-7, FILE- 'UVWIND. DAT ' , STATUS* 'OLD' ) 

CALL  WINDIN ( UWIND ) 

CALL  WINDIN ( VWIND ) 

CLOSE ( UNIT* 7 , STATUS- ' KEEP ' ) 

C 

C  USE  DEFAULT  COMPUTATION  OF  LONG-WAVE  SLOPE: 

SLOPE  -  -1.0 
C 

PHI  -  PI *ANGINC/ 180.0 
DO  20  1-1,21 
DO  10  J-1,25 
C 

C  COMPUTE  WIND  MAGNITUDE  (METERS/SEC)  AND  DIRECTION  (DEG): 

WIND  -  SQRT(  UWIND ( I , J) **2  +  VWIND(I, J) **2  )/100.0 
ANGWD  -  180.0*ATAN2 (  UWIND(I,J)  ,  VWIND(I,J)  )/PI 
C 

C  ANGAZ  DEPENDS  ON  LOOK  ANGLE  AND  WIND  DIRECTION 

ANGAZ  -  ANGLE  -  180.0  -  ANGWD 
C 

C  MAYBE  MUST  ADD  360  DEGREES  TWICE  TO  GET  0  <  ANGAZ  <  360): 

IF  (ANGAZ.LT. 0.0)  ANGAZ  -  360.0  +  ANGAZ 
IF  (ANGAZ.LT. 0.0)  ANGAZ  *  360.0  +  ANGAZ 
C 

C  CONVERT  TO  RADIANS: 

CHI  -  PI* ANGAZ/ 180.0 
C 

C  GET  NRCS  AT  GRID  POINT  (I,J): 

CALL  NRCS (WNM, POL, PHI, CHI, SLOPE, WIND, 
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107  1  I RET, SIGO , SIGDB, CAODB ) 

108  Z ( I / J)  -  SIGO 

109  DZ(I,J)  -  SIGDB 

110  C 

111  10  CONTINUE 

112  20  CONTINUE 

113  C 

114  C 

115  C  Set  text  to  color  index  1 

116  C 

117  CALL  XTCOPT ( ' TCOLOR ',1.0) 

118  C 

119  C  set  up  the  positive,  negative  and  text  colors 

120  C 

121  CALL  CTROPT ( ' PCOLOR ',2.0) 

122  CALL  CTROPT ( 'NCOLOR' ,3.0) 

123  CALL  CTROPT ( 'TCOLOR' ,1.0) 

124  C 

125  C  Set  the  permiter  color 

126  C 

127  CALL  GRIOPT( 'PCOLOR ',5.0) 

128  C 

129  C  ENTRY  EZCNTR  REQUIRES  ONLY  THE  ARRAY  NAME  AND  ITS  DIMENSIONS 

130  C 

131  C  THE  TITLE  FOR  THIS  PLOT  IS 

132  C 

133  C  DEMONSTRATION  PLOT  FOR  EZCNTR  ENTRY  OF  CONREC 

134  C 

135  CALL  WTSTR  (  TX,  TY, 

136  1  '  Normalized  Radar  Cross  Section  (absolute)  ',2,0,0  ) 

137  CALL  EZCNTR  (Z,21,25) 

138  C 

139  CALL  WTSTR  (  TX,  TY, 

140  1  '  Normalized  Radar  Cross  Section  (dB)  ',2,0,0  ) 

141  CALL  EZCNTR  (DZ,21,25) 

142  C 

143  C 

144  C  ENTRY  CONREC  ALLOWS  USER  SPECIFICATION  OF  PLOT  PARAMETERS,  IF  DESIRED 

145  C 

146  C  IN  THIS  EXAMPLE,  THE  LOWEST  CONTOUR  LEVEL  (-4.5),  THE  HIGHEST  CONTOUR 

147  C  LEVEL  (4.5),  AND  THE  INCREMENT  BETWEEN  CONTOUR  LEVELS  (0.3)  ARE 

148  C  SPECIFIED. 

149  C 

150  C  THE  TITLE  FOR  THIS  PLOT  IS 

151  C 

152  C  DEMONSTRATION  PLOT  FOR  CONREC  ENTRY  OF  CONREC 

153  C 

154  C 

155  WRITE(*, 6101) 

156  6101  FORMAT (IX, 'Enter  Area  of  Ocean  Surface  Corresponding  to  ', 

157  1  'Radar  Resolution  (m**2):  ') 

158  READ*,  RDAREA 

159  C 
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160  WRITE (*,6102) 

161  6102  FORMAT (IX, 'Enter  Radar  Cross-Section  of  Target  (m**2):  ') 

162  READ*,  TGTSIG 

163  C 

164  DO  80,  1-1,21 

165  DO  81  J-1,25 

166  AS NR  -  TGTSIG/ (  RDAREA'Z ( I ,  J)  ) 

167  I?  (ASNR.GT.0. 0)  THEN 

168  SNR(I,J)  -  10. 0*ALOG10 (AS NR) 

169  ELSE 

170  SNR(I,J)  -  49.9 

171  WRITE ( * , 6902 )  I,  J,  Z(I,J),  AS NR 

172  ENDIF 

173  81  CONTINUE 

174  80  CONTINUE 

175  6902  FORMAT ( IX ,  ' Z ( ' , 12 ,  ' ,  ' , 12 , ' )  IP,  E14.4, 

176  1  '  AS NR  E14.4,  OP) 

177  C 

178  PRINT  *,  ' - >  Press  "ENTER"  for  SNR  display' 

179  CALL  FRAME 

180  CALL  WTSTR  (  TX  ,  TY, 

181  1  '  Target  SNR  with  Respect  to  Radar  Clutter  ',2,0,0  ) 

182  C-  CALL  CONREC  { Z, 21, 21, 25, -4. 5, 4. 5 , . 3, 0, -1, 0) 

183  CALL  CONREC  (SNR, 21 , 21, 25,-50. 0, +50. 0, 2 . 0, 0, -1 , 0) 

184  CALL  FRAME 

185  C 

186  call  xgend 

187  stop 

188  END 


main  Local  Symbols 


Name  Class  Type  Size  Offset 

SLOPE . local  REAL*4  4  0002 

UWIND . local  REAL*4  2100  0006 

RED . local  REAL*4  24  0172 

GREEN . local  REAL*  4  24  018a 

BLUE . local  REAL*4  24  01a2 

TX . local  REALM  4  Olba 

TY . local  REALM  4  Olbe 

VWIND . local  REALM  2100  083a 

I . local  INTEGER*  4  4  106e 

J . local  INTEGER*  4  4  1072 

TGTSIG . local  REALM  4  1076 

SIG0 . local  REALM  4  107a 

CHI . local  REALM  4  107e 

SLOTS . local  INTEGER*  4  4  1082 

Z . local  REAL*  4  2100  1086 

DZ . local  REAL*  4  2100  18ba 

PHI . local  REALM  4  20ee 

SIGDB . local  REALM  4  20f2 
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main  Local  Symbols 


Name  Class  Type  Size  Offset 

POL . local  INTEGER*4  4  20f6 

ANGLE . local  REAL*  4  4  20fa 

RDAREA . local  REALM  4  20fe 

ANGINC . local  R£AL*4  4  2102 

ANGWD . local  REAL*4  4  2106 

COLS . local  REAL*4  4  210a 

ANGAZ . local  REALM  4  210e 

WIND . local  REALM  4  2112 

WNM . local  REAL*  4  4  2116 

SNR . local  REALM  2100  211a 

AS  NR . local  REALM  4  294e 

IRET . local  INTEGER*  4  4  2952 

NCOLI . local  INTEGER*  4  4  2956 

CA0DB . local  REALM  4  295a 

Symbolic  Constant  Type  Value 

PI . REALM  3 . 1415927E+000 

THPI . REALM  4 . 7123890E+000 

TWOPI  .  REAL *4  6.2831855E+000 

HPI  .  REAL *4  1 . 5707964E+000 

189  SUBROUTINE  WINDIN(WIND) 

190  REAL  WIND(21,25) 

191  CHARACTER* 80  CARD 

192  CHARACTER* 12  HELD 

193  CHARACTER* 4  RFLD 

194  CHARACTER* 1  C(80) 

195  EQUIVALENCE (C, CARD) 

196  EQUIVALENCE (C(  1) , BFLD) 

197  EQUIVALENCE (C< 14 ), RFLD) 

198  C 

199  1  CONTINUE 

200  READ (7 , 7001)  CARD 

201  IF  { BFLD. NE. 'FIELD  TITLE:')  GOTO  1 

202  C 

203  N  -  1 

204  5  CONTINUE 

205  READ ( 7 , 7001 )  CARD 

206  IF  (  (RFLD.NE.'OW  ') .AND. (RFLD. NE. '0E  ')  )  THEN 

207  GOTO  5 

208  ELSE 

209  DO  7  1-1,21 

210  IF  (N.EQ.l)  THEN 

211  READ (7, 7005)  (WIND(I, J) ,  J-  1,  7) 

212  ELSE IF  (N.EQ.2)  THEN 
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213  READ ( 7 , 7005 )  (WIND ( I , J) , J»  8,14) 

214  ELSEIF  (N.EQ.3)  THEN 

215  READ(7 ,7005)  (WIND ( I , J) , J-15 , 21 ) 

216  ELSEIF  (N.EQ.4)  THEN 

217  READ ( 7 , 7005 )  (WIND ( I , J) , J-22 , 25 ) 

219  ENDIF 

219  7  CONTINUE 

220  N  »  N  +  1 

221  IF  (N.LE.4)  GOTO  5 

222  ENDIF 

223  RETURN 

224  7001  FORMAT (A80) 

225  7005  FORMAT ( 6X ,  7F9.2) 

226  END 


WINDIN  Local  Symbols 


Name  Class  Type  Size  Offset 

WIND . param  0006 

I . local  INTEGER* 4  4  295e 

J . local  INTEGER* 4  4  2962 

N . local  INTEGER*  4  4  2966 

CARD . local  CHAR*  80  80  296a 

BFLD . local  CHAR*  12  12  296a 

RFLD . local  CHAR*4  4  2977 

C . local  CHAR*1  80  296a 


Global  Symbols 
Name 

CONREC . 

CTROPT . 

EZCNTR . 

FRAME  . 

GRIOPT . 

NRCS . 

WINDIN . 

WTSTR  . 

XGBEG  . 

XGEND  . 

XTCCOL . 

XTCOPT . 

main . 


Class  Type 

Size 

Offset 

extern  *** 

*** 

*** 

extern  *** 

*** 

*** 

extern  *** 

•  ** 

extern  *** 

*#• 

•  ** 

extern  *** 

*** 

*  *  * 

extern  *** 

*** 

*  *  * 

FSUBRT  *** 

•  *  * 

0736 

extern  *** 

*** 

*** 

extern  *** 

•  *  * 

extern  *** 

*#* 

#*  * 

extern  *** 

*  *  * 

*** 

extern  **• 

*  *  * 

*  ** 

FSUBRT  *** 

*** 

0000 

Code  size  •  097a  (2426) 
Data  size  «  02d7  (727) 
Bss  size  *  29ba  (10682) 


No  errors  detected 
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APPENDIX  D.  Environmental  Sea  Clutter  PACKAGE  SPECIFICATION 


package  Environmental_Sea_Clutter  is 

procedure  NRCS(in  Envir . Surface_Wind  , 

Envir. Swell 

Envir . Precipitation, 

Envir . Salinity  , 

Envir . Air_Temp  , 

Envir . S  ea_Temp  , 

Radar . Frequency  , 

Radar .Polarization  , 

Radar . Inc idence_Angle ; 
out  NRCS  ) ; 

--  Yields  Normalized  Radar  Cross  Section  (NRCS)  of  sea  surface. 

procedure  ERCS(in  Radar . Beam_Specif ication, 

NRCS  ; 

out  ERCS  ) 

--  Yields  radar  cross  section  of  environmental  clutter  for  a  given 
radar.  The  input  information  in  Radar .Beam_Specif ication  can  be 
expressed  in  different  ways,  but  must  be  sufficient  to  determine  a 
"footprint"  over  which  NRCS  should  be  integrated  (e.g.  pulse 
duration  beamwidth,  and  radar  altitude). 

procedure  TSNR(in  Target. RCS, 

ERCS  ; 

out  TSNR  ) ; 

--  Yields  signal- to-noise  ratio  of  a  given  target  with  respect  to 
sea  clutter,  for  the  radar  in  question. 

end  Environmental  Sea  Clutter 
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